Feshbach resonances in ultracold ^'^Li + ^^Na atomic mixtures 



OO 

o 
o 

(N 

1—5 

Oh, 



o 

o 
•I— I 

> 

OS 
(N 

o\ 

O 

o 



X 



Marko Gacesa, Philippe Pellegrini, and Robin Cote 
Department of Physics, U-3046, University of Connecticut, Storrs, CT, 06269-3046 

(Dated: June 6, 2008) 

We report a theoretical study of Feshbach resonances in ®Li + '^'^Na and ^Li + ^"^Na mixtures at 
ultracold temperatures using new accurate interaction potentials in a full quantum coupled-channel 
calculation. Feshbach resonances for Z = in the initial collisional open channel ^Li(/ = 1/2, m/ = 
l/2)-|-^^Na(/ = 1, m/ = 1) are found to agree with previous measurements, leading to precise values 
of the singlet and triplet scattering lengths for the ^'^Li+^^Na pairs. We also predict additional 
Feshbach resonances within experimentally attainable magnetic fields for other collision channels. 

PACS numbers; 



A Feshbach resonance occurs when the kinetic energy 
of two colliding atoms matches the energy of a bound 
level associated with a closed channel Since the 

first experimental observation in a Bose-Einstein conden- 
sate (BEC) of sodium 3, Feshbach resonances have been 
widely used as a tool for controlling interactions in ultra- 
cold atomic gases. One can vary the sign and strength of 
the interactions by tuning the magnetic field near such a 
resonance, allowing the study of ultracold atomic gases 
from the strong interacting to the non-interactingregime 
[|. For instance, BECs of "3Cs d or *^^Rb have 
been realized using this technique to reverse the sign of 
the scattering length from negative to positive. Further- 
more, Feshbach resonances can be used as a spectroscopic 
tool [6.] , and provide a way of making ultracold molecules 
. They allow for the coherent coupling between atomic 
and molecular condensates [§], as well as the creation of 
molecular BECs [9]. More details can be found in Ref. 

m 

Recently, Feshbach resonances have been observed in 
Bose-Fermi mixtures of ^Li and ^^Na |ll| as well as in 



®^Rb and ''^K [13], opening the door to new phenom- 
ena such as boson-mediated Cooper pairing [l3| or the 
formation of heteronuclear molecules Such polar 



molecules may have applications in quantum computa- 
tion [l^ , as well as in the search for the electronic dipole 
moment [l6| . or study of degenerate gases with dipolar 
interactions 171. The most common theoretical treat- 



ment of cold collisions is the coupled-channel calculation 
approach [18i] . although several simplified methods have 
been developed [l^l • These techniques were applied to K- 
Rb mixtures by Simoni et al. [20| and to mixed-isotope 
mixtures of rubidium by Burke et al. 



In this article we present an extensive theoretical study 
of the scattering properties of ^' ''Li and ^'^Na in the pres- 
ence of a magnetic field and at the temperature range 
typical of degenerate gases (i.e. between 300 nK and Ifi 
K). Since Li has both fermionic and bosonic isotopes, 
Li-Na mixtures are of particular interest. We focus our 
efforts on finding the positions and widths of Feshbach 
resonances for various collision channels. This informa- 
tion will be useful for future experiments on cold Li-Na 



mixtures, e.g. to form a BEC of polar molecules. 

For two alkali atoms of relative momentum p and re- 
duced mass fi colliding in a magnetic field, the effective 
two-body Hamiltonian can be written as [22] 
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where Hj is the internal energy of atom j, V 
Coulomb interaction, and the magnetic dipole-dipole 
interaction (often neglected llO|). H^^^ consists of the 
hyperfine and Zeeman contributions, respectively 



H" 
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Here Sj and ij are the electronic and nuclear spin of atom 
j, al^f^ its hyperfine constant (152.1368407 MHz for ^Li, 
401.752 MHz for ^Li and 885.813 MHz for ^SNa), and B 
is the external magnetic field assumed in the z-direction. 

The Coulomb interaction, which depends only on the 
internuclear separation R, can be decomposed into sin- 
glet and triplet contributions 



= Vq{R)P" + Vi{R)P^ , 



(3) 



where Vq (Vi) and P'^ (P^) stands for the singlet (triplet) 
molecular potential and its associated projection oper- 
ator respectively [3]- Molecular potentials were con- 
structed by smoothly connecting the inner part to the 
long-range dispersion expansion 
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where ± are for the triplet and singlet potential respec- 
tively. We used the singlet potential of Fellows [1^ , con- 
structed from accurate spectroscopic data via the inverse 
perturbation approach (IPA), and an ab initio triplet po- 
tential computed using the CIPSI package [2J] and ad- 
justed to match the atomic dissociation energy of Li(2s) 
+ Na(3s) at infinity. For the long range form of both po- 
tentials, we adopted the dispersion coefficients i Cn) and 
exchange energy (Ccx and h) of Fellows (case a) [23|. 
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At ultracold temperatures, s-wave collisions describe 
the scattering process. The collision entrance channel is 
determined by the initial Zeeman states of the atoms. 
A judicious choice of basis is the field dressed molecular 
li]: |a/3), with a = 1,...,6 for ^Li, 
8 for ^Li, and /3 = 1, ... ,8 for ^^Na, all in 
the increasing order of energy. We follow this standard 
labeling throughout the article. In the limit B = 0, it 
reduces to the product of atomic hyperfine states 



hyperfine states 

a = I, 



|/l,mi;/2,r7l2) = |/l,mi)Li «) |/2,m2)Na 



(5) 



where fj 



is the total spin of atom j, and mj its 



projection onto the molecular axis. This basis is suitable 
for description of the collision in the limit of two sepa- 
rated atoms. At smaller separation, the molecular basis 
\SIFMf) becomes more appropriate, with S = si + S2, 
I = ii + 12, F = fi + /2, and Mp its projection. 

Once the bases are defined, one can express Eq.lJT]) in 
matrix form, leading to the following matrix equation 



^F(i?) = ^M(i?)F(i?), 



(6) 



where the coupling matrix M(i?) is defined as 

KJ?^^^^^ = W«,/3)l^1^'(«',/3')), (7) 

and F(i?) is a matrix with columns corresponding to a 
complete set of linearly independent solutions. 

To obtain results with the level of accuracy give n in 
the experimental measurements of Stan et al. we 
needed to include the nuclear Zeeman term "fnijjB and 
the magnetic dipole-dipole interaction term [25| . The 
latter one can be written as 



2 3(R • si)(R • §2) - §1 • §2 
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where ats is the fine structure constant. In the second 
relation, the spin and orbital part are separated, with 
C^{0, 4>) representing a reduced spherical harmonic oper- 

ator and isi(^S2) !j t he second rank tensor operator that 
couples the spins [26|. We took into account the contri- 
bution from by expanding the bases to include the 
total orbital angular momentum I and its projection m;: 
\lmi). We calculated the matrix elements of Eq.® in 
the fully uncoupled basis \lmi; A), where A stands for set 
of quantum numbers {si, m^^ , zi, m^j^ , S2, m^j, i2_^i2}. 
With this basis, (AIF'^IA') takes the simple form (27| 
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Non-zero matrix elements for A/ ~ 0, ±2 and I = I' ^ Q 
couple s and d waves. We included small corrections to 
the energies from higher partial waves up to Z = 4. 

We have also checked that there are no higher partial 
wave resonances {e.g., p-wave, I = I' = 1) for the channel 
1 11). Since the system remains invariant under rotations 
with respect to the axis of magnetic field, the projection 
of the total magnetic quantum number Mp = mi + m2 
is conserved, and it determines the number of coupled 
channels. The Feshbach resonances listed in this work 
are all s-wave resonances (/ = 0). 

We solve Eq.® for different values of magnetic field 
using the multichannel log-derivative method Q to ob- 
tain the S-matrix, from which we extract the radial wave 
function phase shift ri(k) of the initial entrance channel, 
as well as the scattering length a using [2§| 



lim k cot ri(k) 



1 
a 



(10) 



Here, k = ^/2fiE/h is the wave number associated with 
the pair of colliding atoms of relative energy E. 

We first calculated the positions of Feshbach reso- 
nances for the least energetic hyperfine state |a/3) = |11) 
of ^Li-l-^'^Na. There is no inelastic spin relaxation for 
this channel which makes it attractive for trapping since 
the only possible decay channel is through 3-body in- 
teractions. Stan et al. [ll[ have measured the po- 



sitions of three resonances for this state and we used 
their results to calibrate the singlet and triplet poten- 
tials. For each curve, we varied the inner wall by shifting 
the positions of the data points for separations less than 
the equilibrium separation Re according to i?shiftod — 
R+s{R — Re)/{Rc — Re), where s corresponds to the shift 
of the zero-energy classical turning point Re- The best 
agreement with [ll| was obtained with s — 0.06170 a.u. 
(singlet) and s = —0.32878 a.u. (triplet). It may also 
be worth noting that we obtained a very good agreement 
with experimental results only after we included the cou- 
pling with nuclear spin and magnetic dipole-dipole inter- 
action. Neglecting these second order terms resulted in 
agreement up to ±5 Gauss for the measured resonances. 
Our results for the entrance channel |11) and the mag- 
netic field up to 2000 Gauss are shown in Fig. [T] and 
summarized in Table HI Note that four additional Fesh- 
bach resonances were found at higher values of magnetic 
field, at 1097, 1186, 1766 and 1850 Gauss. 

With this set of accurately adjusted potentials, we 
have also determined the singlet (as) and triplet (qt) 
scattering lengths and the energy of the last vibrational 
level (Sj^;^) for both ^Li-l-^SNa and ^Li-f^^^^ mixtures 
(see Table lll| . Stan et al. [11] report a good thermal- 
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FIG. 1: Top: Scattering length for "^Li+^^Na in the entrance 
channel |11). Positions and widths of the resonances are given 
in Table IH Bottom: Zoom on the first three resonances. 
Squares are the experimental positions reported in fH^ . 



TABLE I: ''Li+^^Na - calculated Feshbach resonances (l = 0) 
for the entrance channel |11) and B up to 2000 Gauss. We 
compare with the measured position [ll[, and list the width 
A and the background scattering length Obg. Both Obg and 
A were obtained by fitting to ao = flbg (l — ) . Partial 

waves I = . . .4 were included in the calculation. 



(G) 



pcxp 
^0 



(G) [11] 



A(G) 



flbg (ao) 



746.13 
759.69 
795.61 
1096.68 
1185.70 
1766.13 
1850.13 



746.0 ± 0.4 
759.6 ± 0.2 
795.6 ±0.2 



0.044 
0.310 
2.177 
0.153 
8.726 
0.156 
0.019 



14.003 
13.864 
13.002 
13.902 
12.673 
12.500 
12.499 



ization rate between for the ^Li+^'^Na mixture and ex- 
ploit the effect for efficient sympathetic coohng. Our re- 
sults are in agreement with the estimate based on those 
rates. From the thermalization rate in their experiment, 
the MIT group estimated the ratio of pure triplet elas- 
tic collision between Na-Na (ct^a) and Li-Na (ctab) to 



100 [30|, where aAA — 



''Na-Na 



and 



be aAA/cTAB 

CAB — 47raLi_iMa, with aNa-Na = 85 ao for the triplet 
scattering length of Na, one then estimates |aLi-Na| — 12 
oo for the triplet case. This estimate is also in agree- 
ment with the value of \aT\ ~ 15 oq obtained assum- 
ing a thermalization time Ttot ^ 4t ^ 15 sec, where 
the relationship between r and a is given in [3l|. The 
uncertainty in the scattering length was determined by 
adjusting the inner wall of the potentials to match the 
uncertainty of the resonances (Table [l|. Alternatively, 



TABLE II: '^■'^Li-l-^^Na singlet (S) and triplet (T) scattering 
lengths in units of ao, the last vibrational level, and the cor- 
responding binding energy. See the text for the description of 
the uncertainty. 



Present work 
■^Li+23jg-j^ f^Li+^^Na 



Ref. ^ 

^Li+^^Na 



as 

T 



(MHz) 
(MHz) 



39.7 ±0.5 (1.6) 
36.1 ± 0.3 (4.6) 

47 

12 

-1505 ±3 
-7112 ± 12 



15.9 ±0.3 (1.5) 
12.9 ±0.6 (4.5) 

45 

11 

-1.6 ±0.2 
-5720 ± 16 



39.2 
31.1 



be El,, 

mated value of iJj^gj 



we shifted the triplet potential below the last bound 
state by ±140 MHz [ll|, and then adjusted the inner 
wall of both singlet and triplet curves to match the ex- 
perimental position of the resonances. The last triplet 
bound level was then found with smaller uncertainty to 
5720 ± 16 MHz, as compared to the esti- 
-5550 ± 140 MHz of Stan et al. 
Ill . The same procedure was repeated after varying Cg 
coefficient by ±5%, and the results were within the uncer- 
tainties given in Table [III Note that Feshbach resonances 
are very sensitive to the energy of the last vibrational 
level, which thus can be determined very precisely, while 
the total number of levels can only be estimated since 
it depends on the whole potential, which is usually not 
known accurately enough. For our adjusted triplet state 
of ^Li^'^Na, we found 11 vibrational levels (see Table ITT)) . 

Using our adjusted singlet and triplet potential curves, 
we performed a similar calculation for all entrance chan- 
nels for ^Li+^^Na and ^Li±^^Na. Predicted positions of 
Feshbach resonances are shown schematically in Figures^] 
and [21 It appears that the ^Li±^^Na mixture has sev- 
eral experimentally attainable Feshbach resonances for 
the high-field seeking hyperfine states. Altogether, the 
richness of the scattering properties makes ^'^Li±^'^Na 
mixtures interesting for further exploration. 

In conclusion, we present an extensive study of scatter- 
ing properties in ultracold ^Li+^^Na and ^Li±^'^Na mix- 
tures. A full quantum coupled-channel calculation was 
performed in the field dressed approach, to determine the 
positions and widths of several Feshbach resonances as- 
sociated with different entrance channels. The accuracy 
of our singlet and triplet potentials improved by previous 
experimental measurements allowed us to very accurately 
determine the singlet and triplet scattering length of the 
system, as well as to give an estimated number of vibra- 
tional levels in those potentials together with the energy 
of the last level. 

Z. Pavlovic is gratefully acknowledged for providing 
the original log-derivative program. The authors are 
grateful to O. Dulieu and M. Aymar for providing us 
the LiNa ah initio potentials, and W. Ketterle and V. 
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FIG. 2: Feshbach resonances in ^Li+^'^Na for different coUi- 
sional entrance cliannels. Only channels for which the reso- 
nances were found are shown. Open circles for channel 
were observed experimentally. 
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FIG. 3: Same as Fig. [2] for '^Li+^^Na. 
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